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ABSTRACT 

We have investigated the physical conditions in the narrow-Une region (NLR) of M5 1 using long-sUt spectra 
obtained with the Space Telescope Imaging Spectrograph (STIS) aboard the Hubble Space Telescope (HST) and 
3.6 cm radio continuum observations obtained with the Very Large Array (VLA)^. Emission-line diagnostics 
were employed for nine NLR clouds, which extend 2. "5 (102 pc) from the nucleus, to examine the electron 
density, temperature, and ionization state of the NLR gas. The emission-line ratios are consistent with those 
typically found in Seyfert nuclei and indicate that within the inner near-nuclear region (r < 1") the ionization 
decreases with increasing radius. Upper-limits to the [O III] electron temperature (T^ S 11, 000 K) for the inner 
NLR clouds indicate that photoionization is the dominant ionization mechanism close to the nucleus. The 
emission-line fluxes for most of the NLR clouds can be reproduced reasonably well by simple photoionization 
models using a central power-law continuum source and supersolar nitrogen abundances. Shock-i-precursor 
models, however, provide a better fit to the observed fluxes of an NLR cloud ~2"5 south of the nucleus that 
is identified with the extra-nuclear cloud (XNC). The large [Olll] electron temperature of this cloud (Tg = 
24,000 K) further suggests the presence of shocks. This cloud is straddled by two radio knots and lies near 
the location where a weak radio jet, ~2."5 (102 pc) in extent, connects the near-nuclear radio emission with a 
diffuse lobe structure spanning ^4" (163 pc). It is plausible that this cloud represents the location where the 
radio jet impinges on the disk ISM. 

Subject headings: galaxiss: active — galaxies: individual (M51) — galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 

Photoionization by a central source is widely regarded as 
the dominant process responsible for producing the emission- 
line spectra observed in the Narrow-Line Region (NLR) of 
Seyfert galaxies. Support for these photoionization models 
is provided by their general success in reproducing the ob- 
served ultraviolet (UV) and optical emission line spectra (e.g. 
IFerland & Osterbrock 1986) and by the observed biconical 
morphology present in many Seyfert galaxies, which is sug- 
gestive of collimation by a central ionizing source. 

Many of these Seyfert galaxies also possess radio emis- 
sion, which appears collimated and extended, with sim- 
ilar orientation, on the same scale as the optical emis- 
sion. The radio power has been shown to be correlated 
with the [Olll] luminosity a nd ful l-width half-maximum 
(FWHM) ( de Bruvn & WilsoiJ Il978t |Heckman et al. 198 1 ; 
IWTiittlel[l985t il992i) . High spatial resolution imaging of the 
NLR often illustrates a more detailed correspondence be- 
tween the radio kno ts and th e optical emission- line clouds 
jPogge & de Robertis«il995: .Caoetti et al-1996; f alcke et all 



processes may also contribute. Simple photoionization mod- 
els, in some cases, underpredict the electron temperature in 
the gas (determined from the [Olll] (A4959 + A5007)/A4363 
line ratio) and are not always able to explain strong UV 
emission from intermediate excitation line species, such as 
Cm] and Civ (Ferland & Osterbrock 1986). Thus, the com- 



plex velocity distribution presented by some 


Seyferts iCecil 


1988; Whittle etal. 1988; 


Winse et alJl997; Hutchinss et alJ 


1998) and the aforementioned optical -radio 


correspondence 
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While photoionization models assuming a central ionizing 
source are generally successful in reproducing the observed 
emission-line spectra, there are some indications that other 

' Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy under NASA 
contract NAS5-26555. 

^ The VLA is a facility of the National Radio Astronomy Observatory, 
which is operated by Associated Universities, Inc., under cooperative agree- 
ment with the National Science Foundation. 



suggests that another process, such as shocks, may contribute 
to the observed emission-line spectrum. The role of shocks 
in the NLR provides a key to the understanding the relative 
importance of mechanical versus radiative energy processes 
in the active galactic nuclei (AGN) model. 

We examine the physical conditions influencing the NLR 
kinematics and ionization structure, in particular the relative 
importance of photoionization and shocks within the N LR of 
M5L M51 (NGC 5194) is a nearby (8.4 Mpc, Feldmeie r et alJ 
Il997h . nearly face-on, Sbc grand-design spiral galaxy that 
hosts a weak active nucleus. Early observations attributed 
the ionization structure in the nuclear re gion to p hotoioniza- 
tion b y the active nucleus (.Rose & Searlet 1982: Rose & Cecil 
119831) . However, later observations revealed structures that 
indicated shocks also may be an important source of ioniza- 
tion in this object (Ford et al. 1985; Cecil 1988). In particular, 
shocks were postulated to be the dominant energy mechanism 
exciting the emission-line gas in an extended extra-nuclear 
cloud (XNC) lying ^3" south of the nucleus. 

Several factors led to this conclusion for the XNC. 
iFord et aP Cl985) observed high electron temperatures deter- 
mined from line-ratio diagnostics, emission-line ratios qual- 
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itatively similar to that observed in supernova remnants, a 
limb brightened morphology, and abrupt changes in the spec- 
tral characteristics between the galactic disk and the XNC 
cloud. Further, Cecil ( 1988) found sharply-defined high- 
velocity emission structures that resembled the termini of hy- 
drodynamic flows and the broadest emission lines localized 
along the northern boundary of the XNC. A directed nuclear 
outflow interacting with the disk of M5 1 was postulated as the 
source of the shocks in the XNC. 

Direct evidence for a weak radio jet connecting the nuclear 
radio emission with the southern extra-nuclear emission was 
provided by 4.8 GH z observations obtained with the Very 
Large AiTay (VLA) dCrane & van der Hulsllll992l) . An ex- 
tended lobe structure that is morphologically similar to the 
optical XNC emission also is observed in the radio map at 
this frequency. A lower brightness halo engulfs the nucleus, 
the jet, and regions north of the nucleus. 

Recent Chandra observations (^^1" resolution) by 
l^rashima & Wilson (2001) show an extended X-ray mor- 
phology for the XNC that is similar to both the optical and 
radio data. The X-ray data further support the argument 
for shocks in the XNC. These authors modelled the X-ray 
spectra of this region with a thermal plasma resulting from 
shock-heating due to outflows. Their model yielded a 
X-ray temperature Tx = 6.7 x 10'' K for the XNC, which 
corresponds to a shock velocity of 690 kms"'. This value 
is roughly consistent with the ^5 00 km s~' shock velocity 
previously inferred bv lCecill Jl988h . 

In this paper, we examine the physical conditions influenc- 
ing the NLR kinematics and ionization structure, in particular 
the relative importance of photoionization and shocks within 
the NLR of M51. In §|2]we present the observations, data re- 
duction, and stellar background subtraction. In §|3]we present 
the optical morphology, cloud kinematics, emission-line ra- 
tios, and physical conditions within the NLR gas. In §|4]we 
present our VLA 8.4 GHz radio observations and compare 
them with archival Hubble Space Telescope (HST) Wide Field 
and Planetary Camera (WFPC2) optical images. We discuss 
the photoionization modelling of the NLR clouds in § |5l and 
discuss ionizing shock models in §|6l Our results are summa- 
rized in §0 and our conclusions are presented i n § E| 

Giv en the distance to M51 of 8.4 Mpc aFeldmeier et al.l 
Il997h . 1" corresponds to 40.7 pc. 

2. OBSERVATIONS 

High spatial resolution (0"1, 4.1 pc) long-slit spectra of the 
nuclear region of M51 were obtained with the Space Tele- 
scope Imaging Spectrograph (STIS) aboard HST. The spec- 
tra were acquired through a 52" x 0."2 slit oriented at a sin- 
gle position angle of 166° (E of N). This position angle was 
based on our analysis of an archival //5'77WFPC2 F502N nar- 
rowband image of the near-nuclear region of M5 1 and was 
selected to intersect the nucleus and several of the brightest 
NLR clouds at this wavelength. The low-dispersion grating 
modes G430L (Fig. [0 and G750L (Fig. |3 were used to pro- 
vide continuous wavelength coverage spanning the spectral 
region from 2900 A to 1 fim. The 0"2 wide slit (0."19 as-built 
width) projects onto 3.75 STIS CCD pixels (1.9 resolution el- 
ements) and corresponds to spectral resolutions of 10.3 A for 
G430L and 18.3 A for G750L. The associated velocity resolu- 
tions for this sHt are 827 km s"' at [Oil] AA3727, 615 km s"' 
at [Olll] A5007, and 834 km s"' at Ha. 

To determine the sub-arcsecond velocity structure of the 



emission-line clouds, we observed the spectral region con- 
taining H/3, [Oiii] A4959, and [Oiii] A5007 with the G430M 
(Ac = 4961 A) medium-dispersion grating mode (Fig.|3}- The 
spectral resolution of this data is 1.0 A, which corresponds to 
62kms-i at [Oiii]A5007. 

In addition, we have observed M51 with the VLA in A- 
array at 8.46 GHz (3.6 cm) on 14 August 1999. During the 
data reduction procedure in AIPS, the data were thoroughly 
flagged in order to remove the effects of a severe lightning 
storm that passed over the western arm during the latter part 
of the observations. Four representations of the radio data of 
the central region of M51 are presented in Figure 0] 

A large scale map of the central region (Fig. |4^) displays 
two prominent sources: the nuclear source with an extended 
lobe structure toward the south, and an elongated source at 
27"8 to the northwest at P.A.=335°. The resolution of this 
naturally weighted map is 0"99 x 0."79 and has a rms noise of 
13 /iJy/beam. Higher resolution maps are presented of both 
sources. The naturally weighted map of the nuclear source 
(Fig. I4J5) has a resolution of 0."29 x 0"26 and a rms noise 
of 8 ujy/beam. The uniformly weighted map of the nucleus 
(Fig. 1^) has a resolution of 0"18xO"15 and a rms noise 
of 10 /iJy/beam. While the naturally weighted map shows 
the nuclear source together with the southern jet and the ex- 
tended radio lobe, the uniformly weighted map only shows 
the nuclear source and the inner part of the southern jet. The 
peak flux density of the nuclear source in the two maps is 
210 /j,Jy/beam and 195 /iJy/beam, respectively. 

The discrete northern source in the nuclear region of M5 1 
consists of three compact components and is presented in 
Figure]^. The resolution of this naturally weighted map is 
0."29 X 0"26 and has a rms noise of 8.6 /iJy/beam. 

A log of the optical and radio observations is presented in 
TableU] which includ es the archival H ST/WFPC2 (Ford et J] 
il996t) and VLA 6 cm dCrane & van der Hulsi 19921) data. 

2.1. Data Reduction 

The images and spectra were reduced using the IDL CAL- 
STIS pipeline developed for the STIS Instrument Develop- 
ment Team (IDT). Cosmic rays were removed from the data 
by combining all of the observations in a given spectroscopic 
mode using an iterative sigma-clipping algorithm based on 
a calibrated noise model of the images. After subtracting 
the detector bias and superdark frames, pixel-to-pixel varia- 
tions in the CCD response were removed from the G430L and 
G430M spectra with the standard pipeline flat-field images. 
A contemporaneous flat-field image was used to remove the 
fringing in the G750L spectral image at wavelengths longer 
than 7000 A. 

The data were wavelength caUbrated using exposures of 
an onboard Pt/Cr-Ne spectral line lamp to determine the 
wavelength zero-point offsets to the pipeline dispersion co- 
efficients. Geometric rectification of the images provided a 
constant wavelength along each column of the spectral im- 
ages. Using the pipeline sensitivity tables, the data were flux- 
calibrated in units of erg s"' cm"^ A"' per cross-dispersion 
pixel. 

Residual hot pixels present in featureless regions of the 
post-pipeline processed spectral images were removed using 
an iterative cr-clipping procedure. A 1 x 51 pixel boxcar was 
passed over these regions, and successive filtering iterations 
were performed beginning with 8cr and cleaning down to a 
3(7 threshold. 
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Fig . 1 .— (left) H/3, [O III] A4959, and [O III] A5007 emission lines extracted from the background-subtracted HST/STIS G430L spectral image (see Table[i]for 
details) of the near-nuclear region of M51. The total spatial extent is 8f'6 (350 pc), and the wavelength scale is in the observed frame, (right) Continuum-subtracted 
//57'AVFPC2 F502N naiTowband [OlII] image that has been rotated to the STIS position angle (166°) and resampled to the STIS plate scale (0"051 pixel"'). 
The location of the STIS 52" X 0"2 slit is overplotted on the WFPC2 image. The NLR clouds have been labeled for identification in both images. 



The baclcground present in the medium spectral resolution 
G430M image consists of the near-nuclear stellar continuum 
of the M51 host galaxy, which fills the entire length of the 
STIS slit. Because changes in the continuum shape are small 
(^^5%) over the narrow G430M spectral range (282 A), we 
removed the background using an empirical cross-dispersion 
template of the spectral image in regions without emission 
features. After masking the emission-line clouds, the back- 
ground template was generated by calculating the median of 
each image row. The cross-dispersion template was then sub- 
tracted from each column of the spectral image. 

Nine spatially resolved emission-line clouds were identi- 
fied in the data (see Figures and |2}- Individual cloud spec- 
tra were created by adding adjacent rows along the slit in 
the spatial direction. The number of summed rows (typically 
three pixels, 0" 15) varied from cloud to cloud and was chosen 



based on the spatial extent of the [Oil] AA3727 emission line. 

2.2. G430L and G750L Host Galaxy Subtraction 

Because the active nucleus of M51 is obscured, which is 
expected for a Seyfert 2 galaxy, the AGN continuum is not 
apparent in our lower resolution G430L and G750L spec- 
tra. However, the near-nuclear spectra exhibit a strong stel- 
lar spectrum from the M51 host galaxy. Since the H/3, H7, 
and H(5 spectral features are comprised of both emission and 
absorption lines, accurate measurement and subtraction of the 
stellar components is important for accurately measuring the 
emission-line flux. This is also true for several other emis- 
sion lines (e.g. [Nelll] A3869) that have nearby stellar absorp- 
tion features. Ha absorption is not directly observed due to 
the strong adjacent [Nil] AA6548,6583 emission lines. How- 
ever, the presence of H/3 stellar absorption suggests that Ha 



4 



BRADLEY, KAISER, & BAAN 



STIS G750L 



o i 



SLIT FIDUCIAL 



641 6542 6674 
Wavelength (A) 



STIS G750L 



i 
I 



WFPC2 



F656N 



Cloud 



SLIT FIDUCIAL 



Fig. 2.— (left) [Ol] AA6300,6364, [Nil] AA6548,6583, Ho, and [S II] AA6716,6731 emission lines extracted from the HST/STIS G750L spectral image (see 
Table fn for details) of the near-nuclear region of M51. The total spatial extent is 23'.' 6 (960 pc), and the wavelength scale is in the observed frame. Note the 
extended [Nil] + Ha emission to the north of the nucleus, (middle) Background-subtracted version of the spectral image at left, (right) Continuum-subtracted 
//577WFPC2 F656N narrowband [Nil] + Ha image that has been rotated to the STIS position angle (166°) and resampled to the STIS plate scale (Of'051 pixel"' ). 
The location of the STIS 52" X 0"2 slit is overplotted on the WFPC2 image. The NLR clouds have been labeled for identification in both images. 



absorption also is present. Measurement of the Ha and H/3 
stellar absorption features is especially important because the 
Ha/H/J Balmer decrement is used as our reddening indicator. 
The removal of the stellar background is complicated by the 
presence of variable extinction in the near-nuclear region. The 
variable reddening changes the continuum shape of the spec- 
tra and adds an additional degeneracy to the well known age- 
metallicity degeneracy. If continuum shape alone was used as 
a fit diagnostic, then a highly reddened spectrum could mimic 
an older stellar population. Thus, the stellar population back- 
ground must fit both the continuum shape and the stellar ab- 
sorption lines. 

To determine the stellar background in the near-nuclear re- 
gion of M51, we performed a population synthesis on our 
least extincted (bluest) cloud spectrum in the inner nuclear 
region. This cloud (cloud 6) is 0."43 from the obscured nu- 
cleus. It is less extincted than clouds 4 and 5, and exhibits 
very little differential extinction across its spatial extent. We 
did not expect a stellar population gradient within the near- 
nuclear region (^1", 41 pc) nor did we see evidence in the 



WFPC2 images for localized star formation. Consequently, 
cloud 6 was used to determine our stellar population template 
in the near-nuclear region. To increase the signal-to-noise in 
the spectrum, we combined all four rows (0"2) of cloud 6 
before performing the stellar population synthesis. Assum- 
ing an invariant intrinsic stellar population within the near- 
nuclear region, the only free parameters in fitting the remain- 
ing emission-line clouds are the reddening, E(B-V), and an 
overall multiplicative scale factor 

The first step in the process was to separate the stellar and 
emission-line contributions to the observed NLR spectrum. In 
the case where the emission line occurs in an otherwise fea- 
tureless spectral region, we excluded the emission line from 
the stellar population fit. In regions where the emission from 
the active nucleus is superimposed upon stellar features, we 
simultaneously fit the absorption and emission components 
using the SPECFIT procedure in IRAF^ . These fits were per- 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research in Astron- 
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Fig. 3. — Spectral cuts of the H/3, [Oin]A4959, and [OlII]A5007 emission lines extracted from the background-subtracted HST/STIS G430M spectral image 
(see TablefTlfor details). The clouds are labeled for identification. Clouds with multiple components have the lower-case letters "a" and "b" appended to their 
cloud labels. 



Table 1. M51 Journal of Observations 



Date 


Instrument 


Filter 


Aperture 


Grating 


(A or cm) 


AA 
(A or cm) 


Exposure Time 
(s or hr) 


Optical Data 
















1995 Jan 24" 


WFPC2 


F502N 






5012 


26.9 


1700 


1995 Jan 24= 


WFPC2 


F547M 






5446 


486.6 


860 


1995 Jan 24= 


WFPC2 


F656N 






6564 


21.5 


1800 


1998 Apr 3 


STIS 


Clear 


50CCD 




5852 


4410 


30 


1998 Apr 3 


STIS 


Clear 


52" xO"2 




5852 


4410 


30 


1998 Apr 3 


STIS 


Clear 


52" X 0"2 


G430L 


4300 


2807 


3086 


1998 Apr 3 


STIS 


Clear 


52" x0"2 


G430M 


4961 


282 


2458 


1998 Apr 3 


STIS 


Clear 


52" X 0"2 


G750L 


7751 


4987 


1644 


Radio Data 
















1985 Jan 27'' 


VLA-A 








6 cm 


0.13 cm 


12 hr 


1986 May 9-10'' 


VLA-A 








6 cm 


0.13 cm 


12 hr 


1986 May 11-12'' 


VLA-A 








6 cm 


0.13 cm 


12 hr 


1986 May 15-16'' 


VLA-A 








6 cm 


0.13 cm 


12hr 


1999 Aug 14 


VLA-A 








3.6 cm 


0.04 cm 


10 hr 



"Archival data (Ford et al. 199^^ 

''Archival data tiCrane & van der Hulsi 19921) . 



formed on a line-by-line basis for the spectral region encom- 
passing the absorption and emission lines. A Gaussian profile 
was used for the absorption, while the [O III] A5007 emission- 
line profile was used as a template for the emission. The re- 
sulting absorption fit was inserted into the spectral region of 
these two-component features. 

The spectral synthesis fit wa s then performed on the "stel- 
lar" spectra using a reddened ('Fitzpatrick'11999^ linear com- 
bination of empirical stellar templates (Pickles 1998). The 

omy. Inc., under cooperative agreement with the National Science Founda- 
tion. 



iPicklesI (fT99l stellar li brary was selected over the j Kuruczl 
(?1993i). LSilva&Cornel]l ^992), and Bica & Alloin (1986) li- 
braries due to its complete wavelength coverage, 10 A res- 
olution, extensive stellar types, and the superiority of the 
resulting fit. The 131 spectra in the Pickles stellar library 
spanned all luminosity classes and included both metal-weak 
and metal-rich giants. The stellar synthesis represented the 
best linear combination of the Pickles empirical stellar tem- 
plates plus an overall reddening correction. 

To partially constrain the population synthesis, we used 
the results of J?anagia_et al . (1996J on the dust and stellar 
content in the nuclear region of M51. iPanagia et alJ lfT99^ 
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Fig. 4. — Radio structure in tlie nuclear region of M51 at 3.6 cm witli the VLA-A. (a) (top left) The large scale nuclear structure in a naturally weighted 
map at a resolution of 0"99 X and a rms noise of 13 ^Jy/beam. The contour levels are at [2,4,6,8,10,12,16,20,24,28,32,36,40] X 10 ^Jy/beam. (b) (top 
right) Naturally weighted map of the nuclear region with contour levels at [2,4,6,8,10,12,16,20,24,28,32,36,40] X 7.5 ^Jy/beam. The angular resolution is 
Of'29 X 0"26 and the rms noise is 8 /jjy/beam. (c) (bottom left) Naturally weighted map of the northern source. The contour levels are the same as in (b). The 
peak flux densities of the sources, from north to south, are 64, 68, and 79 /ijy/beam. (d) (bottom right) Uniformly weighted map of the nucleus at a resolution of 
0" 18 X 0" 15. The rms noise is 10 /iJy/beam and the contoiu- levels are the same as in (a). 
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found that the stellar population in the near-nuclear "core" of 
M51 (r < 1") has a main-sequence turnoff around AO (age 
> 400 Myr), while the stellar population in the bulge (r > 1") 
has a main-sequence turnoff near GO (age > 8 Gyrs). Thus, 
for the spectral synthesis on cloud 6 (r = 0"43) we excluded 
O and B main-sequence stars. 

Before performing the stellar population fit, we normalized 
the NLR cloud spectra at 5556 A to match the normaliza- 
tion of the stellar templates. It was unnecessary to broaden 
the absorption features in the stellar templates since the stel- 
lar velocity dispersion in the near-nuclear region of M51 
(62 km s"', Delisle & Hardv 1992) is negligible compared to 
our spectra resolution (615 km s~' at [Olll] A5007) and that of 
the stellar templates (600 km s"' at [Olll] A5007). The NalD 
absorption featu re is conta minated by foregro und Galactic ab- 
sorption ( K ormendv et aL.,1996: Vaceli et al.it 1997) and was 
excluded from the stellar population fit. 

The resulting best fit for cloud 6 (Fig. |5} yields a stellar 
E(B-V) of 0.09 and axl = 0.52 (v = 1506). The lower portion 
of Fig. |5l shows the cloud 6 spectrum after the stellar back- 
ground was removed. The background-subtracted spectrum 
shows very little residual stellar structure with the exception 
of NalD. The NalD feature does not impact our data because 
we do not observe any emission lines in this region of the 
spectrum. 




4000 6000 8000 10000 

Wavelength (A) 

Fig . 5 . — Observed cloud 6 spectrum overplotted with the stellar population 
fit (thick black line) used to subtract the underlying stellar continuum. The 
bottom spectrum shows the results of the background subtraction. 



Comparing the stellar population fit to the near-nuclear 
NLR clouds showed that the stellar composition appears con- 
stant within r < 1" (41 pc, clouds 3-7). Therefore, we fit the 
cloud 6 composite stellar template to the NLR cloud spectra 
within this region by varying the reddening of the stellar tem- 
plate and multiplying it by an overall scale factor 

Because clouds 4 and 5 exhibit variable reddening across 
their spatial extents (5 pixels, 0."25), each of these clouds was 
split into two cloudlets: one that represented the co-addition 
of two rows, and the other the co-addition of three rows. The 
spectra were binned in the spatial direction to both increase 
the signal-to-noise and to perform the measurement on the 
scale of a spatial resolution element (2 pixels, 0."1). All other 
emission-line clouds were summed over their full spatial ex- 
tents (typically three pixels, 0"15) before fitting. 

While the stellar composition appears constant within the 
central radius of (41 pc), the near-nuclear stellar popula- 



tion template provides a poorer fit to the NLR clouds at radii 
> 1". From inspection of the stellar background as a func- 
tion of radius, the relative strengths of the Balmer absorption 
lines decrease while the Call and metal (e.g. TiO and MgH) 
absorption lines increase, indicating an increased contribu- 
tion from an older stellar population outside the near-nuclear 
region. Consequently, we performed a new stellar popula- 
tion synthesis for clouds 2 and 8 (r > 1"). Again using the 
results of Panagia et al. ( 1996), we excluded main-sequence 
stars earlier than GO in the spectral synthesis of clouds 2 and 
8. Figure|6lshows the resulting stellar background fit to cloud 
2. 




4000 6000 8000 10000 

Wavelength (A) 

Fig. 6. — Observed cloud 2 spectrum overplotted with the stellar population 
fit (thick black line) used to subtract the underlying stellar continuum. 



Because we do not observe absorption features in the spec- 
tra of clouds 1 and 9 (r > 2"), we fit the stellar background of 
these clouds with a cubic spline. The background fit to cloud 
1 is shown in Figure^ 




4000 6000 8000 10000 

Wavelength (A) 

Fig. 7. — Observed cloud 1 specfium overplotted with the cubic spline fit 
(thick black line) used to subtract the background continuum. 



2.3. NLR Cloud Velocities 

Cloud velocities and velocity dispersions were measured 
using the bright and unblended [Olll] A5007 fine on the STIS 
G430M spectral image. The central wavelength and FWHM 
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Table 2. Cloud Velocities and Velocity Dispersions 



Uloud 


Radial Velocity 


Velocity Dispersion 




(km s ) 


(km s ) 


1 


-145 (11 


107 (3) 


2 


85 (6) 


202 (18) 


3 


116 (I) 


25 (91 


3a 


217 (61 


138 (13) 


4 


-37 (71 


66 (11) 


4a 


41 (7) 


66 (10) 


4b 


-31 (4) 


331 (13) 


5 


9(1) 


37 (2) 


5a 


-15 (2) 


139 (5) 


5b 


45 (5) 


338 (11) 


6 


30(1) 


22 (3) 


6a 


24 (4) 


221 (13) 


7 


67 (2) 


105 (6) 


8 


53 (7) 


69 (11) 


8a 


-18 (9) 


49 (15) 


9 


34(1) 


74 (3) 



of each [Olll]A5007 line was obtained by fitting a single 
Gaussian plus a linear background to the spectral region ad- 
jacent to the line. Five clouds were not well fit by a single 
Gaussian in the higher spectral resolution G430M data, most 
likely due to a superposition of clouds along the line-of-sight. 
These clouds were fit using multiple Gaussian components. 
We label the additional clouds by appending the lower-case 
letters "a" and "b" to the cloud identifications (see Fig.|3}. 

Several of the emission-line clouds were not centered 
within the STIS slit. We determined the relative position of 
the clouds within the slit by using the HST/WFPC2 F502N 
narrowband image of the near-nuclear region of M5 1 . We ro- 
tated the WFPC2 F502N image to the same position angle as 
our STIS observations (166°) and identified the emission-line 
clouds present our STIS spectra. Gaussian profiles were then 
fit to the NLR clouds in the F502N image to determine the dif- 
ference between the cloud centroid and the center of the slit. 
These pixel differences were converted to velocity shifts and 
were used to correct the cloud radial velocities. At a given po- 
sition along the slit multiple clouds are not spatially resolved 
in the WFPC2 images. Thus, the same slit offset correction 
was applied to clouds with multiple components. 

The intrinsic velocity dispersion of each cloud was cal- 
culated by subtracting, in quadrature, the instrumental line 
broadening of 62 km s"' from the measured FWHM veloc- 
ity widths of the [Olll] A5007 emission Hnes in the G430M 
spectral image. We used the HSTfWFPCl F502N image to 
estimate the instrumental broadening for clouds 2, 3, 5, and 6, 
which do not fully illuminate the slit width. The instrumental 
broadening for these clouds was found to be '--^45 km s"' . 

The G430M [Olll] A5007 cloud velocities and velocity dis- 
persions are presented in Table |2] The uncertainties were 
propagated from the errors in the Gaussian fits and are 
given in parenthesis. A s ystemic velocity of 463 km s~' 
ide Vaucouleurs et alJll99ll) was subtracted from the radial 
velocities to correct for the recessional motion of the M51 
host galaxy. 

2.4. Emission-Line Flux Measurements 

The emission-line fluxes were measured using the IRAF 
routine SPECFIT. Isolated emission lines were measured us- 
ing a Gaussian profile with a linear background. Since we de- 
tected weak residual, broad TiO absorption in the spectral re- 
gion near the [O l] AA6300,6364 lines, a quadratic background 



was fit to the local continuum near these lines and subtracted 
prior to fitting the emission. The [Arm] A7136 flux was mea- 
sured by simple numerical integration across the line profile 
due to the weakness of this emission line and uncertainties 
in the TiO steflar absorption. The G430L [Olll]A5007 pro- 
file was used as a template to deblend the [N ll] + Ha lines 
and the [S ll] AA67 16,6731 doublet in the G750L spectral im- 
ages. Because the G430L and G750L grating modes have 
different spectral resolutions, the width of the G430L [O III] 
profile was allowed to vary in the fits. Using the [Olll] pro- 
file to fit the Ha line allowed for one less free parameter 
and resulted in more reliable fits to the blended [N ll] + Ha 
complex. For cases in which both the Gaussian and [Olll] 
template fit gave reliable results, the Ha fluxes agreed within 
2%. The [Nil] AA6548,6583 doublet Hnes were constrained to 
have the same FWHM velocity widths, fixed relative central 
wavelengths, and an intensity ratio fixed by atomic physics 
(1:2.95). The Ha FWHM was fixed to the value of the 
[Nll]AA6548,6583 lines. Similarly, the [Sll] AA6716,6731 
doublet lines were constrained to have fixed relative central 
wavelengths and identical FWHM velocity widths. 

2.5. Reddening Correction 

We used the observed Balmer decrement of the Ha/H/3 
line ratio to estimate the reddening of each NLR cloud. The 
reddening curve of iFitzpatrickl (11999.) was applied assuming 
an intrinsic Ha/H/3 ratio of 3.1. While case B recombina- 
tion yields Ha/H/3=2.85 for n, = 10^ cm'^ and T, = 10"* K 
( Brocklehurst 1971). an assumed ratio of Ha/H/3=3.1 is more 
appropria te due to colhsi onal excitation of Ha in AGN envi- 
ronments d Ho et al.ll993[ and references therein). The largest 
values of E{B-V) (~0.37, = 1.1 mag) are found for the 
two innermost NLR clouds (4 and 5), which straddle the cen- 
tral dust lane. 

To account for the differential reddening across clouds 4 
and 5, each cloudlet spectrum was separately dereddened be- 
fore summing the cloud flux. The values of E(B-V) and the 
dereddened emission-line fluxes, normalized to H/3, are pre- 
sented in Table |3] The propagated dereddened flux errors are 
given in parentheses. 

3. RESULTS 
3.1. NLR Morphology 

The continuum-subtracted HSTfWFPCl narrowband [O III] 
(Fig.|8t) and [Nil] + Ha images (Fig.lSji) of the near-nuclear 
region of M5 1 suggest the presence of a biconical morphology 
for the inner NLR gas. In the AGN paradigm, ionization cones 
result from the collimation of the nuclear continuum radiation 
field by a dense, optically-thick molecular torus and have been 
pre viously observed in several Seyfert galaxies (Pogge 19881 
lEvansetal.1 1993; Wilson & Tsvetanovlll994 ). Most of the 
emission in the near-nuclear region of M5 1 appears confined 
to an ionization cone with a projected opening angle of '--^74° 
with an axis position angle of ^163° (see Figures|S^ and|SJ)). 

On the southern side of the nucleus, the [N ll] + Ha image 
(Fig.lSji) shows a shell-like structure whose boundary agrees 
well with the boundaries defined by the inner ionization cone. 
The eastern shell boundary is more prominent and better de- 
fined than the western side, which appears to be fragmented. 
Inside of the shell-like structure, we find a hollow cavity that 
is occupied by the southern radio jet. The most probable ori- 
gin of the cavity is that it has been evacuated by the radio 
jet (see § 0}. Thus, the edges of the southern cone structure 



Table 3. Dereddened Emission-Line Fluxes (normalized to H/3) 
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Note. — Fluxes are relative to H/3. Errors are given in parenthesis. 
^'Dereddened H/3 Flux (x 10"'* erg s"' cm-^). 
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may be dominated by shocks expanding laterally from the ra- 
dio jet instead of photoionization from the nuclear continuum 
source. While our single slit position does not intersect these 
boundaries, HST GO program 9147 (Ferruit et al., in prepara- 
tion) will map the emission near the radio jet to determine the 
source of ionization in these structures. 

Additionally, the overall position angle of the optical emis- 
sion is closely aligned with the observe d radio emission at 
3.6 cm and 6 cm (Crane & van der Hulsl fT99l . Recent CO 
imaging of M51 ( S coville et al. 1998) has revealed dense 
molecular gas confined to a disk that is associated with the 
more opaque arm of the X-shaped dust absorption feature that 
crosses the nucleus. This molecular gas is nearly perpendic- 
ular- to the radio jet (Ford et al. 1996; Grillmair et al. 1997) 
and has a density > 10^ cm~^, which is sufficient to collimate 
both the ionizing continuum and the radio-emitting plasma 
((Scoville et al. 1998). This putative molecular torus may also 
provide a reservoir of matter to feed the active nucleus. 

Extensive kinematic modelling by Cecil (1988) indicates 
that the jet axis is oriented ^70° with our line-of-sight. If the 
radio jet and optical bicone axis are similarly collimated and 
nearly co-spatial, then the resulting inclination angle for the 
ionization cone is ^20°, with the southern cone being nearer 

3.2. Kinematics 

We plot the cloud radial velocities and velocity dispersions 
as a function of projected radial distance from the nucleus 
in Figure |9] The amplitude of the NLR cloud radial veloci- 
ties are relatively small and lie in the range from -145 km s"' 
(cloud 1) to 217 km s"' (cloud 3 a). We find the smallest ra- 
dial velocities (|v| < 67 km s"') for the nine clouds north of 
the nucleus (clouds 5-9). Only five of the NLR clouds (1,4, 
4b, 5a, and 8a) have blueshifted radial velocities, while the 
remaining 1 1 clouds are redshifted. The highest velocity dis- 
persions (Av « 340 km s"') are found closest to the nucleus 
(clouds 4b and 5b), and the upper-envelope of velocity disper- 
sions decreases as the projected radial distance increases. 

In Figure|9lwe have overplotted a simple inclined Keplerian 
disk model by assuming the NLR clouds lie in the plane of the 
M51 host galaxy (inclination=20°, eastern side nearer, with 
major axis along RA.= 170°, Tullv 1974) and are undergoing 
gravitational rotation about a central mass of 1 x 10** Mq. It 
is evident that both qualitatively and quantitatively this model 
does not fit the observed radial velocity field of the NLR 
clouds. Most of the clouds north of the nucleus are redshifted, 
however, we would expect them to exhibit blueshifted radial 
velocities if they were located within the galactic disk. The 
lack of any obvious S-shape in the radial velocity distribution 
suggests that the NLR clouds are not undergoing gravitational 
rotation in a circumnuclear disk. 

An alternative explanation for the kinematic structure is that 
the NLR clouds are outflowing (or inflowing) from the nu- 
cleus. Radial outflow motions have been observed in several 
studies of the NLR in Seyfert galaxies (e .g. Kaiser et al. 2000; 
ICrenshaw & Kraemeil200(I iRuiz et all^OOl.) . The relatively 
small observed radial velocities on the northern side of the 
nucleus are consistent with an outflow (or inflow) that lies 
nearly in the plane of the sky. This picture is consistent with 
the clouds lyingnear the front half of the ionization cone (see 
Figures |8t and]^), which we have also inferred to lie nearly 
in the plane of the sky. While our kinematic data cannot con- 
clusively discriminate between outflow and inflow radial mo- 
tions, the presence of the southern radio jet and a possible 
northern counter-jet favor the outflow scenario. 



Clouds 1, 3, and 3a, located on the southern side of the 
nucleus, possess the largest observed radial velocities (|v| > 
116 km s"'). Interestingly, these clouds lie on the same side 
of the nucleus as the radio jet and have nearly coincident ra- 
dio knots (see §|4}. While cloud 1 is consistent with entrain- 
ment by the jet outflow, clouds 2, 3, and 3a are not. For a jet 
orientation of 70° with the line-of-sight (ICecill [T988), these 
clouds should exhibit blueshifted radial velocities if they are 
entrained by the jet. A possible explanation for the larger ra- 
dial velocities of these clouds is that they lie behind the radio 
jet and have offset velocities due to a laterally expanding flow 
away from the radio jet. If there was a weak northern counter 
jet, which has since disrupted, we would expect clouds north 
of the nucleus to be redshifted and those south of the nucleus 
to be blueshifted. We see a very weak signature supporting 
this scenario in both our kinematic and radio observations. 

Assuming that cloud 1 is entrained by the jet outflow 
(70° to the line-of-sight), the resulting deprojected velocity 
of this cloud is 424 kms~'. This velocity is in agreement 
with the previously inferred shock velocity of ^^500 km s"' 
(tCecil 1988) and is roughly consistent with the 690 km s~' 
shock velocity inferred from the Chandra X-ray temperature 
dTerashima & Wilsonl20()l . 

Cloud 1 lies at the northern edge of the XNC near t he lo- 
cation where the radio jet impinges on the radio lobe. Cecij| 
Il988 ) decomposed the XNC into several kinematic compo- 
nents using ground-based (~1" resolution) Fabry-Perot ob- 
servations at 65 km s"' resolution, which is similar to our 
62 kms"' resolution at [Olll]A5007 in G430M. We note 
that his blueshifted XNC component (v = -180 km s"' , Av = 
140 km s"') is in close agreement with our measurements for 
cloud 1 (v = -145kms-i, Av = 104 km s"'). 

3.3. Spectral Classification 

Based upon optical studies of the emission-line ratios, M51 
has been previously classified as a LINER (e.g. Stau ffe] 
1982; Filippenko & Sargent 1985; Carrillo et aljjp^, a 
Seyfert/LINER transition objec t (.Heckman .1980 1). and a 
Seyfert 2 (e.g. iHo et alJ Il997h . Our NLR cloud spec- 
tra exhibit strong emission lines of [O III] AA4959,5007, 
[Nii]AA6548,6583, [Oii]AA3727, and [Sii]AA67 16,6731 
that are typical in Seyfert 2 spectra. 

To cl assify e mission -line objects, 

Baldwin. PhilU ps. & TerlevichI lll981h used a series of 
two-dimensional emission-line ratio diagnostic plots to 
discriminate among various ionization mechanisms. We 
adopted this me thodology and followed the convention of 
IVeilleux & Osterbrock(.1987i. hereafter VO) who chose ratios 
of bright emission Unes that are close together in wavelength 
to minimize the impact of dereddening errors. Figures FTOIll II 
and ini show the three VO diagrams of [O III] A5007/H/3 
versus [Nil] A6583/Ha, [Sll] AA6716,6731/Ha, and 
[Ol] A6300/HQ;, respectively. 

Because the partially-ionized transition zone in nebula pho- 
toionized by young, hot stars is limited in size dStromgrei] 
1939), Hll regions and starburst galaxies typically exhibit 
weaker ionization lines than those observed in LINERs or 
Seyferts. While both LINERs and Seyferts are ionized 
by more energetic sources, LINERs are distinguished from 
Seyferts by their lower excitation, which is indicated in part 
by a lower value of [Olll] A5007/H/3 (Ho et al. 1993). 

In Figures ^| ^\ and El we denote the loci of points 
typically occupied by Hll regions, LINERs, and Seyferts. 
To represent the regions photoionized by OB stars, which 
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Fig. 8.— (a) Continuum-subtracted HSTfWFPCl F502N narrowband [OlII] image and (b) HSTfWFPC2 F656N narrowband [Nil] + Ha image showing the 
biconical emission. The overplotted bicone has a projected opening angle of 74° with an axis position angle of 163° . The location of the 52" X 0"2 slit for our 
STIS observations (P.A.=166°) also is shown. 
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Fig. 9. — Radial velocities (upper panel) and velocity dispersions (lower 
panel) of the NLR clouds as a function of projected radial distance from 
the nucleus. The data points are labeled with the cloud identifications. The 
dashed line represents a simple inclined Keplerian disk model in which the 
clouds are assumed to lie in the plane of the galaxy (i = 20° ) and are under- 
going gravitational rotation about a central mass of 1 X 10** M0. This model 
clearly does not fit the observed NLR kinematics. 



include Hll regions and starburst gala xies, we plotted the 
theoretical Hll model sequence of McCall etalJ (IT985I) for 
T; = 38,500-47,000 K. The M51 emission-hne ratios oc- 
cupy regions consistent with Seyfert galaxies. In general we 
find higher excitation, as measured by the [O lll]/H/3 line ra- 
tios, for the NLR clouds than is typically observed in LINER 
galaxies. Cloud 7, our weakest NLR cloud, has a somewhat 
lower value of [O lll]/H/3, but it lies near the Seyfert/LINER 
boundary. 

Figure [TO] indicates that M51 possesses larger 
[Nil] A6583/HQ! ratios (~ 2.7) than the region typically 
occupied by Seyfert nuclei. We believe this to be an abun- 
dance effect, whereby nitrogen is selectively enhanced in the 




-0.5 0.0 
log ([Nil] X6583/Ha X6563) 

Fig. 10.— [Oiii] A5007/H/3 vs. [Nil] A6583/Ha diagnostic plot for the 
emission-line clouds. The loci of points typically occupied by Seyferts and 
LINERs ai'e denoted by the colored boxes. The green hne represents the 
theoretical Hll region track of McCal l et alj IMsh . The arrow indicates the 
direction of the reddening correction. 



near-nuclear region of M51 by factors of ^ 3-4.5 (see § l5.2t . 

3.4. Ionization State of the Gas 

The ionization state of the NLR gas can be exam- 
ined through the use of the [Olll] A5007/H^ ratio and 
the abundance insensitive [Olll] A5007/[Oll] AA3727 
and [Olll]A5007/[Ol]A6300 Une ratios. The 
[Olll] A5007/[Oll] AA3727 line ratio is roughly propor- 
tional to the ionization parameter for radiation-bounded 
clouds and is relati vely independent o f the shape of the 
ionizing continuum ( Penston et a 1119901) . 

Figures O-c show [Olll] A5007/H^, 

[O III] A5007/[O II] A A3727, and [O iii] A5007/[O i] A6300, 
respectively, as a function of projected radial distance from 
the nucleus. In general, these line ratios are larger near 
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Fig. 11.— As in Fig. [To| but for [OlII] A5007/H/3 vs. 
[Sll] AA6716,6731/Ha. 
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Fig. 12.— As in Fig.[To| but for [OlII] A5007/H/3 vs. [Ol] A6300/Ho. 



the nucleus and decrease with increasing radius out to 
'^l". As a comparison, we have also overplotted three r~" 
(n = 0.5,0.6, 1.0) curves, which qualitatively agree with the 
data, on Figures I13h-c. The declining line ratios indicate 
that the excitation of the NLR gas decreases with increasing 
radius within the near-nuclear region. As an exception, cloud 
1 possesses the largest value of [O III] A5007/H/3 (even larger 
than the innermost clouds 4 and 5), which may indicate the 
presence of an additional source of ionization (e.g. shocks) in 
this cloud. 

3.5. Electron Density 

In the low-density regime, the [S ll] A A67 16,6731 doublet 
is a well-known electron density diagnostic that is relatively 
insensitive to temperature (Osterbrock 1989). This diagnostic 
can be applied in regions where the electron density is less 
than the critical density of [S ll] (3.3 x 10^ cm"-'). Above this 
density the lines become collisionally de-excited. Unfortu- 
nately, the [S II] AA67 16,6731 lines are severely blended in 
our data. The wavelength separation of these lines (14.4 A) 
is less than the spectral resolution of the G750L grating mode 
(18.3 A). We deblended the [S ll] AA6716,6731 doublet us- 
ing the [O III] A5007 line as a template, however, the result- 



ing flux errors are large (Fig. 11311). We note that for a [S ll] 
A6716/A6731 ratio near unity, a 10% error in each of the 
emission lines results in nearly a factor of two uncertainty in 
the electron density. 

To calculate the electron density of each NLR cloud (Ta- 
ble EJ, we used the "temden" routine of the IRAF NEBU- 
LAR package, w hich uses a five-level atom approximation 
(HE aw & Dufoudll995 ). For clouds 2-9, we assumed an elec- 
tron temperature Te = 10"^ K, which is typical for photoion- 
ized gas in the NLR of Seyfert galaxies. For cloud 1, we used 
an electron temperature of 2.4 x 10^ K, which was c alculated 
from the observed [Olll] emission lines (see § I3.6> . The re- 
sulting electron densities lie in the range between 490 cm""* 
and 1150 cm"^, however, no trend in the gas density as a 
function of radius can be inferred given the large flux er- 
rors. The average [S ll] electron density within the NLR gas 
is 770 cm"-'. Cloud 1, located at the largest radial distance 
(2"5), possesses the lowest electron density of n^ = 490 cm"-'. 

We compare our results with those of Rose & Searle ( 198:^ 
who adopted an average value of [Sll] A6716/A6731 = 1.0 
for the entire nuclear region of M51 (61. "4 x 2" slit). These 
authors assumed = 8 x 10^ K and calculated an electron 
density of 630 cm"^. If we adopt this electron temperature, 
we find an average electron density of 690 cm"^, which is in 
very good agreement with their value. Further, our results 
are similar to those of Rose & Cecil ( 1983) who found in a 
5" circular- aperture an average [Sll] A6716/A6731 = 0.93, 
corresponding to an electron density of 880 cm"^ for T,, = 
7620 K. Our calcula ted electrq r i den sities are also consistent 
with those found by Ford et aLlJl98 5) who obtained electron 
densities in the range 640-720 cm"^ for the nucleus of M5 1 
and 370-490 cm"^ for the XNC. 

3.6. Electron Temperature 

The electron temperature can be estimated from the ob- 
served [Oiii] (A4959 + A5007)/A4363 ratio. The difficulty 
in using this diagnostic is that the [Olll] A4363 emission line 
generally is very weak in Seyfert spectra. For our NLR cloud 
spectra, we observed [Olll]A4363 only in cloud 1. Using 
the NEBULAR package, we calculated an electron temper- 
ature of 24,000 K for this cloud (assuming an electron den- 
sity of 490 cm""*). Electron temperatures of this scale typi- 
cally are not observed in the O zone of photoionized gas 
(Kraemer et al. 1998), providing further evidence that shocks 
contribute to the ionization structure of this cloud. 

Because the [O III] A4363 line was too weak to be detected 
in the other NLR clouds, we derived upper limits for the elec- 
tron temperature of the other clouds by numerically integrat- 
ing the flux over the strongest noise feature in the spectrum 
near 4363 A (Table 0}. The resulting electron temperature is 
£ 1 1,000 K for the three innermost clouds (4-6) and cloud 8. 
These temperatures effectively rule out collisional ionization 
as the dominant source of excitation in these clouds as elec- 
tron temperatures < 20, 000 K a re too low for this mecha- 
nism to be important (Wilson*1 993^ . 

Clouds 3 and 9 have electron temperatures < 17,500 K, 
which also suggests that the ionization of these clouds is pri- 
marily due to photoionization. The upper limits for the tem- 
perature of clouds 7 and 2 are 19,000 K and 21,800 K, respec- 
tively. It is most likely that these represent weak constraints 
and the large upper limits are due to the lower signal-to-noise 
of these clouds, which are fainter than the inner NLR clouds 
and cloud 1 . While it remains that the larger upper-limit to the 
electron temperatures could possibly indicate a greater role 
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Fig. 13.— (a) [OlII]A5007/H/3 (b) [OlII] A5007/[OlI] AA3727 (c) [OlII] A5007/[Ol] A6300 and (d) [SlI] A6716/[SlI] A6731 for the nine NLR emission-line 
clouds plotted as a function of radial distance from the nucleus. 



Table 4. Physical Conditions in the NLR Clouds Along P. A. 166° 



Cloud 


Radius" 


lie 




L(H/3) 


Ionized Mass 




(",pc) 


(cm-3) 


(K) 


(lO'*^ ergs-') 


(M0) 


1 


2.53 (103) 


490 


24,000 


2.7 


50 


2 


1.07 (44) 


510 


< 21,800 


3.0 


60 


3 


0.83 (34) 


680 


< 17,500 


6.5 


90 


4 


0.18 (7) 


850 


< 10,200 


31.1 


350 


5 


-0.16 (7) 


620 


< 9,900 


30.2 


470 


6 


-0.43 (18) 


840 


< 11,000 


9.8 


110 


7 


-0.75 (31) 


1090 


< 19,200 


2.0 


20 


8 


-1.05 (43) 


1150 


< 11,100 


3.7 


30 


9 


-1.91 (78) 


670 


< 15,800 


1.7 


30 


Total 








90.7 


1210 



'Positive values denote clouds south of the nucleus. Negative values denote clouds north of the nucleus. 



for shocks, the MAPPINGS II shock models (see §|6} provide 
a poorer fit to these clouds than our photoionization models. 

3.7. Ionized Mass 

The mass of ionized gas in each NLR cloud can be esti- 
mated from the observed electron density and luminosity of 
the H/3 recombination line. For a completely ionized gas with 
solar helium abundances, the total ionized mass is given as 
jOsterbrockil98l 

where mu is the mass of a hydrogen atom, L(H/3) is the 
observed luminosity (erg s"') of the H/3 recombination line, 
a^^l is the Case B effective H/3 recombination coefficient 



(3 X lO-i"* cm^ s-i at 10"* K, 'Osterbrock"198^, n, is the 
electron density, and hvnp is the energy of an H/5 photon. 
We derived the H/3 luminosities by correcting the observed 
H/3 emission-line fluxes for the distance to M51 of 8.4 Mpc 
(Feldmeier et al. 1997). The electron densities were taken as 
those derived in § l3.5l for the [S ll] emitting region. The result- 
ing ionized masses of each NLR cloud are shown in Table |3 
The total ionized mass of the NLR clouds in our 52" x 0."2 slit 
is 1 .2 X 10^ Mq, with most of this mass located in the three in- 
nermost and brightest H/3 clouds. For comparison. Ford et a^ 
Jl985h derived a total ionized mass of 4.4 x 10^ Mq within a 
larger 2" x 2" region centered on the nucleus. While our slit 
position was chosen to intersect the bright [O ill] clouds strad- 
dling the nucleus and the bright NLR cloud located near the 
jet terminus, additional structures, such as two bright clouds 
north of nucleus and faint extended circumnuclear emission. 
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are included in Ford's measurement. 

3.8. Ionizing Flux 

The observed H/3 luminosity can also be used to calculate 
the to tal number of recombinations per second as JPetersonI 
fT997l) 

Q{H) = ^"^^^^ = 2.11 X 10'2l(H/3) photons s"' , (2) 

where as is the Case B hydroge n reco mbination coefficient 
(2.6 X 10"'^ cm^ s-i at 10"^ K, Osterbrock 1989). 

For the nine NLR clouds located within our 52" x 0."2 slit 
along position angle 166°, the total H/3 luminosity is 9.1 x 
lO-'^ erg s"' (see Table|4}- Using the above equation, we de- 
rive a hydrogen recombination rate of Q{H) = 1 .9 x 10''° pho- 
tons s~' . This number provides an estimate of the total num- 
ber of hydrogen ionizing photons produced per second by the 
active nucleus by assuming ionization equilibrium, whereby 
the total number of hydrogen ionizing photons produced by 
the central source equals the total number of recombinations. 

Because the observed H/3 luminosity was summed over 
our single slit position, we have not accounted for all of 
the NLR H/3 emission. As such, Q{H) represents a lower- 
limit to the rate of ionizing photons produced by the central 
source. Within a 2" x 2" aperture, which as discussed previ- 
ously contains emission structures not included in our narrow 
slit. Ford et alJ (tl985) reported a nuclear H/3 luminosity of 
2.8 X 10^^ erg s~' . This results in a total rate of hydrogen ion- 
izing photons of 5.9 x lO^'* photons s~'. These calculations 
show the relative weakness of the activ e nucleus of M51. For 
comparison, Schulz & Komossa ( 1993) derived a total rate of 
hydrogen ionizing photons of 3 x 10^'* photons s~' for the cen- 
tral source of the luminous Seyfert 1 galaxy NGC 4151. 

4. RADIO OBSERVATIONS 

The existence of a radio jet and radio lobe on the southern 
side of the nucleus of M51 wa s first detected using VLA ob- 
servations at 4.8 GHz jCrane & vanderHulst 1992). These 
authors found a sinuous radio jet, extending ^^3", that con- 
nected the weak nuclear radio source with a diffuse radio 
lobe, which has b een identified with the optical XNC feature 
jFbrd et al.ll985 l. The radio lobe was dominated by a bright 
arcuate, ridge-like structure extending approximately 5". 

Our VLA-A 8.4 GHz observations (see Fig.|4} confirm the 
presence of the previously observed radio structures in the 
near-nuclear region of M5 1 . The radio maps show a weak, 
unresolved radio core with a flux density of 330 /iJy/beam. 
The southern radio jet spans 2"3 (94 pc) and decreases in 
luminosity along its spatial extent, reaching a minimum of 
~2-4 /iJy/beam. The position angle at the base of the jet is 
158°. The radio jet terminates near a diffuse, extended lobe 
structure, however, it is not connected to the radio lobe at the 
lowest contour level. It appears that the radio jet connects to 
the southern radio lobe near the location of cloud 1, however, 
because of the low signal-to-noise of our radio observations 
we are unable to determine the precise location of the im- 
pingement. Although our res olution at 3.6 cm is be tter than 
the earlier 6 cm data (Crane & van der Hulslll 19921) . a com- 
parison of the peak flux in the jet at these two wavelengths 
suggests that the spectral index is approximately 3.1. The to- 
tal east- west extent of the southern radio lobe is 5. "3 (220 pc). 

The nuclear source of M5 1 is also surrounded by a 2" ex- 
tended radio structure that is elongated along a position angle 



of 169°. To the north-northwest of the nucleus, this struc- 
ture turns into a ridge that extends approximately l."5 and 
could signify a bent northern counter-jet. Our 3.6 cm data 
do not exhibit the northern loop structure observed at 20 cm 
by Ford et al. ( 1985) and at 6 cm by C rane & van der HulstI 
(fl992). Because this feature is only weakly detected at 6 cm, 
it is likely not present in our maps due to the lower signal- 
to-noise of our radio data and is probably not a result of the 
frequency difference. 

In addition to the nuclear source structures, we find an elon- 
gated triple radio source at 27 ."8 (1130 pc) north-northwest 
of the nucleus at a position angle of 335°. At the distance 
of M51, the spatial length of this structure is 86 parsec. The 
position angle of this northern source is within 14° of the posi- 
tion angle of the nuclear elongation, and it is possible that this 
triple structure is a remnant of past jet activity. However, it 
could also result from a weak radio galaxy in the background. 

The centerline of the radio jet and its brightness distribution 
reveal a substantial non-linearity that is suggestive of bends 
and knots in the jet. We identify eleven radio knots along the 
entire north-south structure. The location and the flux den- 
sity of these radio knots are presented in Table |5] The earlier 
4.8 GHz radio maps only displayed two distinct hotspots. 

Our VLA 8.4 GHz data reveal a slightly different jet orien- 
tation as compared to the earlier low-resolution 4.8 GHz data 
(Crane & van der Hulst 1992). Specifically, the position an- 
gle of the southern radio jet in our 8.4 GHz maps (observed 
in 1999) lies ~10° to the east of the centerline of the ear- 
Her4.8 GHz jet structure (observed in 1985/1986). Addition- 
ally, the brightness distribution of the jet and southern radio 
lobe structure appears different than the earlier observations. 
In particular, the southern lobe structure is rather patchy and 
does not exhibit the clear arcuate, emission-ridge structure 
that was observed in the earlier 4.8 GHz data. 

A potentially intriguing explanation for the discrepancy in 
the radio structures in the south would be oscillatory struc- 
ture in the jet. The total radio power in the nuclear region 
of M51 (2 X 10'^ W Hz"') is at the very b ottom of the dis- 
tribution of active galaxies (ICondonlll992^ . and the relative 
weakness of the radio jet makes it susceptible to instabili- 
ties that can cause the jet to precess. Oscillatory structure 
may be the result of precession at the nucleus itself or from 
Kelvin-Helmholtz (KH) instabilities resulting from the veloc- 
ity differences between the gas inside and outside the jet. The 
KH instability results in helical structure for the jet that shows 
an exponentially growing oscillation with distance along the 
jet. However, the deduced shift in position angle of a pre- 
cessing jet is relatively large and would result in a precession 
period that is too short. The mis-orientation of the low- and 
high-resolution jets could be due to a limb-brightening effect, 
but this would require rather unrealistic frequency-dependent 
opacity variations across the width of the jet. Alternatively, 
the apparent shift could simply result from the signal-to-noise 
limitations of our data and t he uniformly weighted map of 
ICrane & van der Hulst* ("19921). 

In Figures nn and [T5l we present flie HSTfWFPCl [Olll] 
and [N ll] + Ha narrowband images superposed with the VLA 
8.4 GHz radio continuum contours. In order to align the radio 
core with the assumed location of the obscured nucleus, the 
radio image (0"29 x 0."26 resolution) was shifted relative to 
the optical images (0."09 resolution) by Aa = 0"15 and A6 = 
-1''06. The location of the nine optical emission-line clouds, 
the eleven radio knots, and the sht for the STIS spectra also 
have been indicated. 
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Table 5. Radio-Optical Correlations 



Radio Cloud 


Radius" 
(", pc) 


Radio Flux 
(/.tjy/beam) 


Optical Cloud 


Radius" 
(", pc) 


H/3 Flux 
(10""' erg s"' cm"^) 


Comments 


JNortn 














KiU 


-l.S (/i) 


z.o 








Adjacent to Cloud 9 


R9 


— 1 2 (49"! 


2.2 




—1 05 f43^ 


4.43 


AHiarpnt to ClmiH 8 


R8 


-0 s ns"! 


1.6 


Cloud 7 


-0 75 f31^ 


2.32 


Adjacent to Cloud 8 


R7 


—0 6 (24) 


2.8 






1 1 62 


AHiii^f^nt tc\ Chimin r\ 
/^UltlCClll LO \^10LIU U 


R6 


-0.3 (12) 


3.6 


Cloud 5 


-0.16 (7) 


35.73 


Adjacent to Clouds 5 & 6 


Nucleus 





155 











South 














R5 


0.6 (24) 


3.6 


Cloud 3 


0.83 (34) 


7.66 


Cloud 3 at edge ofR5 


R4 


1.0 (41) 


1.6 


Cloud 2 


1.07 (44) 


3.54 


See morphology note'' 


R3 


1.4 (57) 


2.6 








See morphology note*" 


R2 


1.7 (69) 


2.6 








See morphology note'' 


Rl 


2.3 (94) 


3.1 


Cloud 1 


2.53 (103) 


3.16 


RO & Rl straddle Cloud 1 


RO 


2.7 (110) 


2.6 


Cloud 1 


2.53 (103) 


3.16 


RO & Rl straddle Cloud 1 



"Positive values denote clouds south of the nucleus. Negative values denote clouds north of the nucleus. 
''Morphologically, the region between R4 and Rl appears as though it may have been evacuated by the radio jet. 




53.1' 53.0' 52.9' 52.8' 52.7' 52.6' 52.5' 1 3" 29"' 52.4' 
Right Ascension (J2000) 



Fig. 14. — Continuum subtracted //57yWFPC2 F502N narrowband [OlII] image superimposed with our VLA 8.4 GHz radio contours. The contour levels 
are [2,3,4,5,6,7,8,9,10, 11,12, 16,20,24,28,32,36,40] X 5.0 /jjy/beam. The location of the optical NLR clouds and radio knots have been indicated as well as the 
position of the STIS 52" X 0"2 slit. 



The radio structures are elongated with a position angle 



close to that of the optical NLR emission observed in the HST 
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53.1= 53.0= 52.9= 52.8= 52.7= 52.6= 52.5= 1 S*" 29"' 52.4' 
Right Ascension (J2000) 



Fig. 15. — Continuum subtracted HSTIWP¥C2 F656N narrowband [Nil] + Hq image superimposed with our VLA 8.4 GHz radio contours. The contour levels 
are [2,3,4,5,6,7,8,9,10,11,12,16,20,24,28,32,36,40] X 5.0 /jjy/beam. The location of the optical NLR clouds and radio knots have been indicated as well as the 
position of the STIS 52" X Q'.'l slit. 



images. The STIS 52" x 0."2 slit falls on most of the inner 
part of the elongated nuclear radio structure. Six of the nine 
optical NLR clouds located within our slit lie near the pro- 
jected position of radio features. At the nucleus, the extended 
[Nil] + Ha and [Olll] emission falls within the outer radio 
contours. The overall morphological agreement of the opti- 
cal and radio emission at 8.4 GHz is very suggestive. We find 
that the optical emission-line clouds, in general, lie adjacent to 
the peaks of the radio emission. This spatial anti-coincidence 
suggests that the radio jet and its structural components do af- 
fect the optical emission regions and are associated as an en- 
ergy source with the optical emission-line clouds. However, 
there are a couple of radio knots within the radio jet/optical 
bicone region that do not have nearby optical knots and thus, 
not all of the radio knots are anti-correlated with an emission- 
line cloud. While the radio resolution is 3 times worse than 
the optical image, it appears that the radio and optical knot 
anti-coincidence would still be present in a higher resolution 
radio map. South of the nucleus, it appears that the radio jet 
may have evacuated a conical region (see Fig.ll5>. In Table|5] 
we provide a more detailed comparison of the locations and 



properties of the radio knots with the H/3 emission from the 
optical clouds. 

The [Olll] and [Nil] + Ha //primages show that cloud 1, 
located 2."5 south of the nucleus, falls within the radio lobe at 
the northern edge and coincides with a clear radio hotspot at 
location RO-Rl . This hotspot could be interpreted as the entry 
point of the jet into the radio lobe. However, earlier observa- 
tions by Cecil ( 1988) associate higher velocity dispersion gas 
east of this point with the entry point of the jet. Higher signal- 
to-noise radio observations are necessary to confirm the pre- 
cise location at which the radio jet impinges the radio lobe. 

The WFPC2 F656N [N ll] + Ha image shows an extended 
shell-like structure that surrounds the region where the ra- 
dio jet is currently passing through. This shell is particularly 
prominent on the eastern side of the source, while it is frag- 
mented on the western side. It is intriguing to interpret this 
shell structure as the boundary of a cavity region that has 
been evacuated by the radio jet. The excitation of this shell 
structure could be due to radiation emerging from the energy 
dissipation within the jet. Alternatively, the excitation of the 
shell could result from lateral expansion shocks generated by 
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the twisting jet. 

5. PHOTOIONIZATION MODELS 

We have generated photoionization models to determine if 
the NLR emission-line spectra are consistent with photoion- 
ization by a non-thermal central continuum source. The pho- 
toionization models were co mputed using versio n 94 of the 
photoionization code Cloudy (Ferl and et alJl 998'^). The mod- 
els were parameterized by the shape of the ionizing contin- 
uum, the dimensionless ionization parameter U, and the total 
hydrogen density hh- The dimensionless ionization parameter 
U is defined as 

1 Ly 

A-TTr^nnc J^^ hv 

where r is the distance from the inner (illuminated) face 
of the cloud to the central continuum source, hh is the to- 
tal hydrogen density, Iivq is the ionization energy of hydro- 
gen (13.6 eV), and L^, is the frequency-dependent luminos- 
ity of the nuclear ionizing continuum. The models assume 
a constant density cloud with a plane-parallel geometry and 
emission-line photon escape from the illuminated face of the 
cloud. 

5.1. Model Input Parameters 

Photoionization models with solar abundances cannot re- 
produce the relatively strong observed [N l] AA5200 and 
[Nil] AA6548,6583 niti-ogen lines in M51. We demonstrate 
this graphically in Fig. where we have plotted the pho- 
toionization model grids for both solar and four times solar 
nitrogen (4 N0). The best-fitting models require selective 
enhancement of nitrogen in the range of 3-4.5 times solar. 
The numerical abundances relative to hydrogen for a model 
with a factor of three overabundance of nitrogen are: He = 
1.0 X 10"', C = 3.55 X lO-"*, N = 2.8 X lO"'', O = 7.4 x 
Ne = 1.17x10-4, S = 1.62xl0-^ and Ar = 3.98x10"^ 
JOrevesse & Anders 19891 
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-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 

log ([Nil] >.6583/Ha J.6563) 

Fig. 16.— [Oni]A5007/H/3vs. [No] A6583/Ha for the NLR clouds over- 
plotted with our Cloudy photoionization model grids. The solar abundance 
model grid is shown in magenta and is labeled with log U = —2.5 to -3.5 and 
log uh = 2 to 4 cm"^. The red grid shows the photoionization models with a 
4 X enhancement of the nitrogen abundance. 

We also generated photoionization models that included the 
presence of dust grains within the line-emitting region. How- 
ever, when compared to our observed optical emission lines. 



these models are very similar to models without dust grains. 
Ultraviolet spectra and/or emission lines involving elements 
commonly found in dust grains can be useful for constrain- 
ing models with dust. The lUE spectrum (10" x 20" elliptical 
aperture) of the nuclear region in M5 1 only shows the pres- 
ence of a blended Sim] A 1892 and C m] A 1 909 emission line 
attributed to Hll regions (Ellis et alJ ll982h . We investigated 
whether the absence of ClvAA1549 from the WE spectrum 
could be used to constrain dust models. For our highest ex- 
citation photoionization model, the ClvAA1549 reddening- 
corrected line flux 1 x 10"'^ erg s"' cm"-) would be un- 
detectable in the M51 lUE spectrum. Therefore without a 
more sensitive UV spectrum, we are unable to constrain mod- 
els with dust grains. While the photoionization models that 
follow do not include the presence of dust grains within the 
NLR line-emitting region, we cannot rule out such models. 

Using published M51 continuum flux measurements, we 
constructed a broadband spectral energy distribution (SED) to 
define the shape of the ionizing continuum. We assumed the 
SED was a power-law of the form /v oc where the spec- 
tral index, a, was defined piecewise for different frequency in- 
tervals. Numerous photoionization models were generated to 
investigate the effect of changing the SED slope in all regions 
of the spectrum, in particular the UV to X-ray wavelength 
regime. As one would expect, the models were most sensitive 
to the spectral index in the extreme ultraviolet (EUV) region 
of the spectrum. We started with aEuv = 1 -29 (hi/ > 13.6 eV), 
which was derived from the published data, and constructed 
several photoionization model grids by varying aEuv over the 
range 1-2. We found the best-fitting photoionization mod- 
els occurred for aEuv = 1-35. In addition, we tested a model 
with the presence of an ultraviolet blue-bump in the ioniz- 
ing continuum. However, this model produced a poorer fit 
to our observed line fluxes, especially for [N l] AA5200 and 
[Oi]AA6300,6364. 

The photoionization models were generated using Cloudy's 
optimize command, where we allowed the ionization param- 
eter, U, and total hydrogen density, rtfj, to vary for each 
cloud. As previously discussed, our early photoionization 
model grids without optimization indicated that enhanced ni- 
trogen abundances were required to fit the observed emission- 
line spectra. Therefore, the relative nitrogen abundance also 
was allowed to vary in the modelling procedure. 

5.2. Model Results 

The photoionization model results are presented in Fig- 
ureEl where we plot the ratio of the model predictions to the 
observed emission-line fluxes. With the exception of cloud 1, 
the vast majority of the emission lines are reproduced reason- 
ably well (within a factor of two) by the models. Given the 
simplicity of the models and the limited number of input pa- 
rameters, this agreement suggests that photoionization is the 
dominant ionization mechanism in most of the NLR clouds. 

High-excitation [Nev]A3346 and [Nev]A3426 lines were 
detected in the spectra of the two NLR clouds closest to 
the nucleus (clouds 4 and 5). The highly-ionized [Nev] 
emission lines have a large ionization potential (97.1 eV) 
and were not produced by our single-component, optically- 
thick (radiation bounded) photoionization models. The pres- 
ence of these high-ionization emission lines indicates that an 
additional optically-thin (matter bounded) component is re- 
quired in these clouds, similar to that found by studies of 

other Seyferts (,Binette et al 1996: Komossa & Schulz. 19971 

lKraemeretalJl2000h . The two-component models were gen- 
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erated by assuming that all of the [Nev] emission arises in 
the high-ionization matter-bounded component and that lower 
ionization lines (e.g. [Olll], [Oil], [Nil], and [Sll]) are pro- 
duced in the radiation-bounded component. As before, these 
models were optimized in Cloudy by allowing the ionization 
parameter, total hydrogen density, and nitrogen abundance to 
vary. For the matter-bounded components, we also varied the 
total hydrogen column density, which controls the depth of 
the cloud. The resulting two-component models for clouds 4 
and 5 are shown in Figure Single-component, radiation- 
bounded models were used for all of the other NLR clouds. 
The model parameters used to generate the jphotoionization 
models for each cloud are presented in Table|6l 

The model results show the presence of a gradient in both 
the ionization parameter and the nitrogen abundance. Exclud- 
ing cloud 1, which is poorly fit by the photoionization mod- 
els, the ionization parameter is larger for the clouds closest 
to the nucleus and decreases with increasing radial distance 
out to cloud 9 ('^2", 80 pc). This result confirms the gen- 
eral trend observed earlier in the [Olll] A5007/[Oll] AA3727 
and [O III] A5007/H/? emission-line ratios and is consistent 
with the radial dilution of a central ionizing radiation field in 
which the gas density decreases more slowly than r"^. Inter- 
estingly, the photoionization model results also show a gra- 
dient in nitrogen abundance over the NLR region (Fig. I18> . 
The nitrogen abundance appears to be larger closer to the nu- 
cleus (^4.5 Nq) and decreases with increasing radius, reach- 
ing ^3 N© for the outer clouds. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Radial Distance from Nucleus (arcseconds) 

Fig. 18. — Nitrogen abundance (relative to solar) of the optimized Cloudy 
photoionization models for the NLR clouds. 



Peimbert ( 1968) was the first to claim that solar abundances 
could not explain the observed emission-line ratios in the 
near-nuclear region of M5 1 . His work suggested that nitrogen 
was overabundant in the nucleus of M51 (r < 3."5) by a fac- 
tor ofjibout five relative to solar. Later emission-line studies 
hvlRlTse & Searle (1982) andlFord eTalJ ll1985[ -19 Nq for 
the XNC) also inferred an enhancement of nitrogen in the cir- 
cumnuclear region ofM51. Ourlai-ge [Nil] AA6548,6583/Ha 
ratios (2.8 - 6.2) confirm that nitrogen is enhanced in the NLR 
ofMSl. 

Nitrogen overabundances on the scale inferred by us 
in M51 have also been claimed in studi es of man\ 
Seyfert galaxies ( Storchi-Beramann & Pasto rizalll989t ll^ 
iDoDita & Sutherland 1996; Storchi-Bergmann et al.' Il99g 
LINERS (Binette 1985; FiHcDenko & Terlevich 1992|, and 
starburst galaxies (Cozioletal. 1999). In particular, 
JStorchi-Bergmann & Pastoriza 1990) found that a range of 
nitrogen abundances, from 0.5 to 5 times solar, provided a 



better fit to the the NLR emission-line ratios from a sam- 
ple of 180 Seyfert 2s and LINERs. Of these galaxies, ap- 
proximately one-third presented a nitro gen overabundance in 
the nucleus ( Storchi-Bergmannl 1199 Ih . These results sug- 
gest that nitrogen is overabundant in active nuclei when- 
ever [Nil] AA6548,6583/Hq; > 1.5. The enrichment of ni- 
trogen found in the nuclei of Seyferts may indicate that 
the N/O abundance gradient found in the disks of spiral 
galax ies (|Searle 1971; Smith 1975; Vila Costas & Edmundi 
Il993h continues into the nuclear regions. Generally speak- 
ing, the selective enhancement of nitrogen may be caused 
by secondary CNO nucleosynthesis, in which nitrogen form s 
from preexisting carbon and oxygen d Thurston et al.llT996l) . 
and may be distributed by stellar winds from massive stars 



and may be distributed by s 
jCid Fernandes et al.ll992l) . 



5.3. Model Discrepancies 

While the overall agreement of the simple photoionization 
model predictions and the observed optical spectrum gener- 
ally is good, there are a few discrepancies. The poorest match 
occurs for the [Oil] AA3727 emission line, which is generally 
overpredicted by the models. For three of the NLR clouds 
this discrepancy is within a factor of 2.5. The largest mis- 
match (a factor of —3) occurs for cloud 6. Although it is 
possible that our measured [Oil] fluxes are low as a result 
of dereddening errors, additional dereddening would change 
the emission-line fluxes for the entire spectrum. In particu- 
lar, the [Nelll] AA3869,3968 lines, which are at a comparable 
wavelength to [O ll] AA3727, would also increase in strength 
counter to the model predictions. Given the generally good 
agreement of the models with the dereddened line fluxes, the 
discrepancy between the [Oil] AA3727 model prediction and 
the dereddened [O ll] AA3727 flux does not appear to be a re- 
sult of dereddening errors. 

One method to reduce the predicted [O ll] AA3727 
emission-line fluxes is to increase the total hydrogen den- 
sity, hh- [Oil] AA3727 has a relatively low critical density of 
3.3 X 10^ cm"^, above which the lines are collisionally sup- 
pressed. To examine this effect, we have generated photoion- 
ization models with increased total hydrogen density. While 
the [On] AA3727 lines are indeed suppressed, the higher den- 
sity models provide a poorer overall fit to the observed spec- 
trum. In particular, [N l] and [S ll] are suppressed and many 
Unes, including [Nelll] A3869, [Olll] A5007, [Ol] A6300, and 
[S III] A9532, are overpredicted for these models. It is pos- 
sible that our simple assumption of constant density clouds 
is invalid and that more complex models involving density 
stratification may provide a better fit for [Oil] A A372 7. Or 
perhaps, the nitrogen overabundance described in § l5.2l is cor- 
related with an oxygen underabundance. 

We also find modest discrepancies for the [S III] and [S ll] 
emission lines. [S III] is slightly overproduced (by a factor 
of —1.4) and [Sll] is slightly underproduced (by a factor of 
— 1 .4) by the models. A general problem with photoionization 
models is that they overpredict the strength of the [S III] A9069 
and fSlll]A9532 Unes (Ho et al. 1993; Kraemeret al. 2000). 
lOsterbrock et alJ J 19921) suggested that this discrepancy may 
be improved with underabundances of sulfur, however, we 
note that t he dielectronic recombination rates for sulfur are 
uncertain (l AUet al.ll!99TI) . While the [S III] lines are over- 
predicted, the [Sll]A6716 and [Sll]A6731 lines are a some- 
what underpredicted. This in turn may suggest the pres- 
en ce of a overabundance of sulfur, whic h has been detected 
bv lStorchi-Bergmann & Pastorizal (119900 in a sample of 180 
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Table 6. Photoionization Model Parameters 





log U 


Hydrogen Density 


Column Density 

/ 1 n20 \ 

(10 cm ) 


JNitrogen Abundance 


Note 






(10^ cm"^) 


(Nq) 




1 


-2.83 


4.6 


6.1 


3.4 


RB 


2 


-3.27 


2.3 


1.9 


3.0 


RB 


3 


-3.07 


0.7 


2.6 


4.0 


RB 


4a 


-3.13 


2.6 


2.8 


4.3 


RB 




-1.52 


1.0 


0.8 


4.3 


MB 


5" 


-2.90 


1.8 


4.6 


4.3 


RB 




-1.08 


0.2 


1.8 


4.3 


MB 


6 


-3.26 


6.3 


2.5 


4.3 


RB 


7 


-3.39 


4.8 


1.7 


3.6 


RB 


8 


-3.21 


2.5 


2.2 


3.9 


RB 


9 


-3.39 


5.1 


1.8 


3.1 


RB 



Note. — RB: radiation bounded; MB: matter bounded. 
"90% RB + 10% MB. 
^$0% RB + 20% MB. 



Seyfert and LINER galaxies. Clearly, altering the sulfur abun- 
dances will not improve the fit for both [S III] and [S ll]. How- 
ever, the discrepancy in these lines is relatively modest. 

As noted earlier, the photoionization models provide a poor 
match to the observed spectrum of cloud 1 . This cloud also 
possesses the largest [Olll]/H/3 ratio, indicating the presence 
of high ionization. The proximity of this cloud to the entry 
point of the radio jet into the lobe suggests that shock excita- 
tion may be an important source of ionization. Previous stud- 
ies (Ford et al. 1985: Cecil 1988: Terashima & Wilson 2001) 
have provided considerable evidence that shocks contribute to 
the ionization structure of the XNC in M5 1 . 

6. SHOCK IONIZATION MODELS 

While it is widely held that the dominant source of energy 
input into the NLR gas is photoionization by a nonstellar nu- 
clear continuum (Koski 1978; Osterbrock 1989), collisional 
shock wave heating may provide an additional so urce of ion- 
ization in the NLR (Morse et al. 1996: Bicknell etalJll998t 
IWilson & R avmond 1999). Shock waves generated by su- 
personic radio ejecta would compress and heat the ambient 
gas and may contribute, a t least locally, a sig nificant fraction 
of the ionizing photons ( iDopita & Sutherland .1995. .1996t 
iBicknell et al. 1998). These EUV photons are generated in the 
hot {T ^ 10^ K) postshock cooling region and diffuse both up- 
stream and downstream. The photons diffusing downstream 
influence the ionization of the shock recombination region, 
while those diffusing upstream ionize the preshock medium. 

Dopita & Sutherland ( 1995, Il996h have created a low- 
density grid of fast radiative shock models using the code 
MAPPINGS II. The models assume a steady-state one- 
dimensional radiative flow with roughly solar abundances. 
iDopita & Sutherland! (1 1 9951 [l996) emphasized the impor- 
tance of the EUV photons that diffuse upstream into the 
preshock gas and produce a high ionization precursor region. 
We refer to the MAPPINGS II shock models that include 
the postshock cooling and photoionized gas as "shock-only" 
models, while those that include contributions from both the 
shocked gas and the ionized precursor region as "shock + 
precursor" models. Current, publicly available, shock mod- 
els are limited because they represent a one-dimensional 
model of a more coniplex t hree-dimensional shock structure 
jDopita & Sudierlandll996l) . 

The two most important parameters in controlling the emer- 



gent shock spectrum are the shock velocity, Vj, and the mag- 
netic parameter, B/n'/^. The shock velocity controls the 
shape of the ionizing continuum, and the magnetic param- 
eter controls the effective ionization parameter o f the post- 
shock recombination zone (IDopita & Sutherlandll995..1996l) . 
The MAPPINGS II model grids are computed for shock ve- 
locities in the range of v, = 150-500 km s"' and magnetic 
parameters B/n^l^ = 0-4 ^Gcvcv'l^. Magnetic parameters 
B/n^l^ ~ 2-4 /iG cm^/^ correspond to the equipartition of 
the magnetic and thermal pressures in the preshock medium 
f)opita & Sutherland 1996). 

We have compared our observed NLR line ratios with the 
MAPPINGS II shock models to determine if shocks play a 
significant role in the excitation of the NLR gas. As demon- 
strated earlier with the Cloudy photoionization models, there 
appears to be an overabundance of nitrogen in the near- 
nuclear region of M5 1 . Similarly, we find that the solar abun- 
dance MAPPINGS II model grids cannot reproduce the large 
observed [Nll]/Ha ratios o f the NLR clouds in M51. For 
a sample of active galaxies. iDopita & Sutherland! ( 1199^ ob- 
served [N Il]/Ha ratios that were ~0.4 dex larger than the so- 
lar abundance MAPPINGS II shock model grids. These au- 
thors likewise suggested that this discrepancy was due to an 
enhancement of the nitrogen abundance and generated new 
shock models with a N/O abundance of two times solar. They 
found a nearly one-to-one increase in the strength of the [N ll] 
lines with the scaled abundance. Thus to approximate the ef- 
fect of nitrogen enhancement in M51, we scaled the MAP- 
PINGS II nitrogen emission-Une fluxes ([Nl] and [Nil]) by 
the relative nitrogen overabundance. 

The MAPPINGS II shock-only models are ruled out for 
all of our NLR clouds because they fail to predict the ob- 
served [Oii]AA3727/[Oiii]A5007 and [Oiii] A5007/H/3 line 
ratios. The M51 [Oil] AA3727/[Olll] A5007 line ratios lie 
in the range 0.1-1.1, while the shock-only models predict 
[On] AA3727/[Olll] A5007 > 3. This is not a surprising result 
as the shock-only models generally lie in the region occupied 
by LI NER galaxies ( Dopita & Sutherland^! 995; Dopita et al] 

However, the MAPPINGS II shock-nprecursor 
models can reproduce the observed 

[Oll]AA3727/[Olll]A5007 and [O III] A5007/H/3 line ra- 
tios for shock velocities in the range 300 - 500 km s"' 
(Jx ~ 1.3-3.5 X lO*" K), consistent with results for several 
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other Se yfert ga laxies dPopita & Sutherlan d! 1995HAllen et all 
11299; Ferruit et al. 1999b' Evans et al. 1999). Unfortunately, 
for shock velocities in this range there is considerable overlap 
of the shock+precursor and photoionization model grids for 
most of the optical emission lines. As such, optical emission 
lines generally are unable to provide conclusive discrim- 
ination between shock an d photoionization models (e.g, 
i Dopita & Sutherland! fT995t lAUen et alJ fT999t lEvans et alJ 
Il999h . However, shocks provide a rich, colhsionally-excited 
UV spectrum, many of whose lines are predicted to be 
much stron ger than t hat fo und by photoionization models 
JDopita & ^S utherlan dl 1 19961 [Allen et alj Il998h . Because 
of this, lAllen et al.i lll998h emphasized the importance of 
ultraviolet emission lines, such as NlllA991, ClvAA1549, 
He II A 1640, Cm] A 1909, and Cll]A2326, in distinguishing 
between shock and photoionization models. Unfortunately, 
the ultraviolet emission in M5 1 is very weak, and a sensitive 
UV spectrum of the near-nuclear region of M51 does not 
exist. Furthermore, the available lUE UV spectrum only 
shows a blended Sim] A 1892 and Cm] A 1909 emission line, 
which has been attributed to Hll r egions located within the 
large 10" x 20" aperture ( Ellis et alJlHl). 

Without UV emission lines, the best diagnos- 
tic we have available to differentiate between the 
shockH-precursor and photoionization models is 
[Olll]A5007/H/? vs. [Ol]A6300/Ha (Fig. [TSi. As shown 
in Fig. shockH-precursor models with shock velocities 
Vj > 300 km s"' predict a stronger [Ol] A6300/Hq! ratio 
than our photoionization models. Overall, our observed 
[O l] A6300/Ha line ratios appear to be more consistent with 
the photoionization model grids. Excluding clouds 1, 9, and 
7, the shockH-precursor models overpredict the strength of the 
[O l] A6300/Ha ratio by factors of 2-6, effectively ruling out 
the shock+precursor models for these clouds. 
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Fig. 19.— [Oiii]A5007/H/3vs. [O I] A6300/Ha for the NLR clouds over- 
plotted with the MAPPINGS II shock-only (green), shock+precursor (blue) 
and Cloudy model grids with 4x solar nitrogen (red). The MAPPINGS II 
model grids ai'e labeled with shock velocity in units of km s"' (shock-only: 
150—500 km s"', shock+precursor: 200-500 km s"') and magnetic parame- 
ter B/n'''^ (0-4) in units of /iG cm^^^^. The Cloudy photoionization model 
grid spans log U in the range -3.5 to -2.5 (bottom to top) and log hh from 2 
to 4 cm"'' (left to right). 

While the shock+precursor models are able to predict the 
strength of the [O l] A6300/Hq; ratio in clouds 7 and 9, they 
fail to provide a good fit for all of the observed emission 
lines. For cloud 7, the best overall fit shock+precursor 



model is for = 250 km s"' and B/n'/^ = 1 ^Gcm^l^. 
However, in this model the strength of the high ionization 
He II A4686/H/3 ratio is underpredicted by a factor of five. 
The He II A4686/H/5 prediction is improved by increasing 
the shock velocity to 350 kms"', but the [Ol] A6300/H/3 
and [Nl] AA5200/H/5 ratios are overproduced by a factor 
of -3 and the [O ill] A5007/H/3 and [O ii] AA3727/H/3 ra- 
tios are overproduced by factors > 2. Turning to cloud 9, 
the best overall shock+precursor model fit is also for = 
250 kms"' and B/n'/^ = 1 /iG cm^/l While this model 
provides the lowest x^, many of the emission lines, includ- 
ing [Neiii]A3869, [Oiii]A5007, and [Nii]A6583 (even with 
3.1 N0), are underproduced by factors of ^2, and Hell A4686 
is underproduced by a factor of four. The fit appears to be 
influenced by the relatively strong observed [0 11] AA3727, 
[Nelll]A3869, [Nl]AA5200, and [Ol]A6300 emission lines. 
Excluding these lines from the calculation, we find the 
best shock+precursor model parameters are = 350 km s"' 
and B/n^/^ = 1 /iG cm^/^. While this model provides an im- 
proved fit for most of the emission lines, [Nl] AA5200 is over- 
produced by a factor of four and [Oil] AA3727, [Ol] A6300, 
and [Nelll]A3869 are overproduced by ^2. There does not 
appear to be a particular combination of and B/n^l^ that 
provides a good fit to all of the emission lines in cloud 9. 
While shocks cannot necessarily be ruled out in clouds 7 and 
9 based on our optical emission-line data, we find that our 
photoionization models provide a better overall fit to these 
NLR clouds than the shock+precursor models. 

6.1. Cloud 1 

In contrast, we find that shocks are important for cloud 1. 
As discussed in § 15.31 photoionization models do not provide 
a good fit to observed optical spectrum of cloud 1 (see Fig ll7t . 
However, the MAPPINGS II shock+precursor models pro- 
vide a reasonable fit to the observed emission lines in this 
cloud. The best fit shock+precursor models have parameters 
in the range = 400-450 km s"' andB/n'/^ = 2-4^G cm^/^ 
with scaled nitrogen abundances of '^3.5. In Fig.|20|we plot 
the ratio of the v, = 450 km s"' and B/n'/^ = 4 ^G cm^/^ 
shock+precursor model predictions to the observed emission- 
line fluxes in cloud 1. For comparison, we have also over- 
plotted the optimized Cloudy photoionization model results 
for this cloud. In particular, the shock+precursor models 
produce a better match to the [Nl] AA5200, [Ol] A6300, and 
[Ol] A6364 lines than the photoionization models, which un- 
derproduce these emission lines by factors of ^3. In Ta- 
ble we compare the optical spectrum of cloud 1 with the 
shock+precursor model predictions for the full range of shock 
velocities and magnetic parameters. 

For a fully ionized gas with solar helium abundances, the 
predicted postshock temperature can be estimated as 7^ = 
1.4 X 10^ 

^jioo where Vjioo is the shock velocity in units 
of 100 km s-i jHoUenbac h & McKe3 ll979l) . For shock ve- 
locities in the range Vj = 400 - 450 km s~' , we derive post- 
shock temperatures in the range from 2.2 - 2.8 x 10'' K. These 
values are roughly consistent with the X-ray temperature of 
Tx = 6.7 X 10'' K obtained from Chandra observations of the 
XNC iTerashima & Wilson 2001), providing further evidence 
that shocks are important for the excitation of cloud 1 , which 
lies near the northern edge of the XNC. 

We do not find evidence of a broad [O III] A5007 velocity 
dispersion for cloud 1 . However, because this line is produced 
primarily in the shock precursor region, this may be expected. 
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Table 7. Line Ratios for Cloud 1 and MAPPINGS II Shock+Precursor 
Models 
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4.44-4.90 
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0.56-0.80 


0.87-1.33 
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1.17-2.86 


1.28-2.95 


1.34-2.94 


Ha A6563 


3.09(1.51) 


2.99-3.02 


2.92-2.95 


2.90-2.91 


2.90-2.91 


2.89-2.92 


[Nil] A6583 


14.33 (6.85) 


2.01-4.31 


2.40-6.66 


3.44-8.41 


3.77-8.69 


3.93-8.66 


[Sll]A6716 


1.90 (0.92) 


0.29-0.56 


0.58-0.91 


0.81-1.21 


0.78-1.16 


0.82-1.11 


[Sll]A6731 


1.71 (0.84) 


0.25-0.41 


0.49-0.70 


0.58-1.14 


0.57-1.23 


0.59-1.27 


[Arm] A7 136 


0.49 (0.24) 


0.08-0.11 


0.11-0.20 


0.10-0.23 


0.09-0.23 


0.09-0.22 


[S III] A9069 


0.96 (0.44) 


0.24-0.33 


0.25-0.49 


0.24-0.56 


0.24-0.60 


0.25-0.65 


[S III] A9532 


1.47 (0.64) 


0.59-0.80 


0.60-1.18 


0.59-1.37 


0.59-1.45 


0.60-1.58 



Note. — The range of values in columns 3—7 correspond to the range of magnetic parameters B/n}l^ - 0—4 /xG cm-'^^. The nitrogen lines have been scaled 
by 3.5 to approximate the effect of an overabundance of nitrogen. 
"Normalized to H/3 unless otherwise indicated. 
''Dereddened fluxes with errors shown in parenthesis. 
^ [O III] A4363/[0 III] A5007 electron temperature ratio. 
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Fig. 20. — Ratios of the MAPPINGS II shock+precursor model predictions 
and the optimized Cloudy photoionization model predictions to the dered- 
dened emission-line fluxes of cloud 1 . The MAPPINGS II shock-^precursor 
model is for i', = 450 km s"' and B/n^l^ = 4 /iG cm^^^ with the nitrogen lines 
scaled by 3.5 to approximate the effect of nitrogen overabundance. 



For lines which are produced in the postshock region, such as 
[On] AA3727, [Olll] A4363, and [Nev] A3426, we only have 
lower spectral resolution data and are therefore unable to de- 
termine if the velocity widths are broader for these lines. 

The largest discrepancies between the MAPPINGS II 
shockH-precursor model and the observed spectrum of cloud 
1 are the underprediction of the [S ll] and [O III] A4363 emis- 
sion lines. As the [S III] A9069 emission line also is slightly 
underpredicted, it appears the model fit could be improved 
with a modest overabundance (^^ 2x solar) of sulfur An 
overabundance of sulfur previously has been found in a sam- 
ple of Seyferts and LI NERs and may be correlated with the 
nitrogen enhancement iStorchi-Bergmann & Pastorizalll989l 



119901 ISchmitt et all 1 19941) . However, the discrepancy with 
the [S II] lines may be a result of the incomplete dielectronic 
recombination rates for sulfur ( Ali et al.-1991 1 

While the [Olll] A4363/[Olll] A5007 ratio is under- 
predicted in cloud 1, this appears to be a general 
feature of the MAPPI NGS II shock-nprecursor models 
iPopita & Suther land 1995). This is especially true for shock 
velocities Vj > 300 kms~', in which the majority of the 
[Olll] emission is produced in the lower-temperature pho- 
toionized precursor rather than the hot postshock region 
( Sutherland etal. 1993; Dopita & Sutherland 1995, 1996). 
The [Olll] A4363/[Olll] A5007 prediction is improved with 
shock-only or low-velocity shock-nprecursor models, how- 
ever, these models provide a much worse overall fit to the 
ensemble of emission lines (see Table 0. In any case, the 
MAPPINGS II shockH-precursor models match the observed 
optical spectrum of cloud 1 better than the photoionization 
models (see Fig.l20>. 

6.2. H[3 Luminosity 

Correcting for the distance to M51 (8.4 Mpc), the observed 
cloud 1 H/? luminosity is 2.7 x 10^^ erg s~'. The H/9 lumi- 
nosit y predicted from the sh ock-Hprecursor models is given by 
iT^ita & Sutherll^lfT995h as 

, 2.41 



L(H/3) = A no 



7.44 X 10" 



100 km s" 



9.85 X 10" 



2.28' 



(4) 



lOOkms" 

where A is the shock area (cm"^) and no (cm"^) is the preshock 
density. Because the cloud 1 shock is vie wed nearly edge on 
(see § 13. U . we follow the examples of Ferruit et al. (1999al) 
and Schmitt et al.. C2002.) and estimate the shock area by as- 
suming a line-of-sight depth equal to the slit width (0."19, 
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8 pc). For this cylindrical symmetry we derive a shock area 
of 4.5 X 10^^ cm"^ (47 pc^). We estimate the preshock den- 
sity using the shock compression factors for the [Sll] zone 
computed by Ferruit et al. ( 1999a). For the best-fitting shock 
velocities in the range from 400-450 km s"' these compres- 
sion factors are 46.4-50. Given the calculated [Sll] electron 
density of 490 cm"^ for cloud 1, the resulting preshock densi- 
ties are in the range from 9.8- 10.6 cm""'. For shock veloci- 
ties Vj = 400-450 km s"' , the shock-nprecursor model predicts 
H(3 luminosities in the range L{H(3) = 2. 1 - 2.6 x lO"'* erg s"' , 
which is approximately 78-96% of the observed H/3 luminos- 
ity of cloud 1 . This suggests that, within the uncertainties in 
the preshock density and the estimated shock velocity, shocks 
are the dominant excitation mechanism in cloud 1 . 

7. SUMMARY 

We have examined the physical conditions in the NLR 
of M51 using archival HSTfWFPCl narrowband images, 
high spatial resolution, long-slit HST/STIS spectra, and VLA 
8.4 GHz radio continuum maps. Our velocity measurements 
of the NLR clouds show that the radial velocity distribution 
is inconsistent with gravitational rotation in a circumnuclear 
disk. The kinematics are consistent with radial outflows from 
the nucleus that lie nearly in the plane of the sky. On the 
southern side of the nucleus, the radial velocities of the NLR 
clouds appear to be influenced by the presence of the southern 
radio jet. 

Through the use of standard line-ratio diagnostics, we have 
shown that the emission-line ratios in the NLR clouds are in 
agreement with those typically found in Seyfert galaxies. The 
ionization state of the gas, as probed by the [Olll]/H/3 and 
[Olll]/[Oll] line ratios, decreases radially from the nucleus 
out to r ^ 2". This is consistent with photoionization from 
a nonstellar nuclear continuum source in which the density 
falls off more gradually than r~^. The ionization is larger for 
cloud 1 (r = 2"5) and appears to be caused by the presence of 
shocks. 

The [Sll] A6716/A6731 flux ratios indicate that the aver- 
age electron density within the NLR gas is ^770 cm"^. Given 
the large error bars for the [S ll] fluxes, no radial trend can be 
inferred for the electron density. Cloud 1 is the only cloud 
in which we detected [Olll] A4363. The derived [Olll] elec- 
tron temperature is 24,000 K. This, coupled with our radio 
observations, suggests that shocks are an important heating 
mechanism for this cloud. Upper limits to the [O III] electron 
temperature for clouds 4, 5, 6, and 8 (7^ < 11,000 K) indi- 
cate that collisional heating is less important for these clouds. 
These electron temperatures are consistent with photoioniza- 
tion. 

The optical and 3.6 cm radio emission both exhibit an ex- 
tended, knotty emission structure along a position angle of 
170°. Our 3.6 cm data are in agreem ent with the earlier 
6 cm observations of lCrane & van derHulst(.,1992) . Both ob- 
servations exhibit a weak southern jet, ~2f'5 in extent, that 
connects with a diffuse lobe structure south of the nucleus. 
The well-known XNC optical structure south of the nucleus 
appears to be coincident with the southern radio lobe. The 
higher spatial resolution 3.6 cm data exhibit several knots 
of emission, indicating that there is structure in the jet. In 
general, these radio knots lie adjacent to the observed op- 
tical clouds, suggesting that the radio jet may contribute to 
the kinematics or ionization structure of the optical emission. 
Cloud 1, at the northern edge of the XNC, appears to lie near, 
but west of, the entry point of the jet into the southern ra- 



dio lo be. Recent Chandra observations (iTerashima & Wilsonl 
also detect X-ray emission coincident with the XNC. 

At the assumed entry point of the radio jet into the diffuse 
lobe structure, we find both radio and optical emission. The 
optical emission, cloud 1, is straddled by two radio knots, RO 
and Rl. The proximity of cloud 1 with the jet-lobe entry point 
and radio knots RO and Rl suggests that colUsional excitation 
may be important for this cloud. 

Cloud-by-cloud photoionization models were generated to 
determine the relative importance of photoionization for each 
cloud. With the exception of cloud 1, the photoionization 
models provide good fits to the NLR emission-line clouds. 
However, our models require supersolar nitrogen abundances 
to match the dereddened [Nl]AA5200, [Nll]A6548, and 
[N II] A6583 line ratios. This result is supported by other stud- 
ies of Seyfert galaxies that also required selective enhance- 
ment of nitrogen by factors up to five times solar to fit the 
observed nitrogen emission. In addition, we find a strong 
gradient in the nitrogen abundance in the near-nuclear re- 
gion. The nitrogen abundance is ^^4.5 Nq for the near-nuclear 
clouds (r ^ 0."2) and decreases for increasing radial distance 
to ^3 N© for the outer clouds (r ^ 2"). 

We also find a gradient in ionization parameter, which 
agrees with the result inferred from the observed emission- 
line ratios. We detect the highly-ionized [Nev]A3346 and 
[Nev]A3426 emission lines in clouds 4 and 5 (r ^ 0"2), 
which indicates the presence of a high-ionization matter- 
bounded component. We successfully modelled these clouds 
with two-component models, which included contributions 
from both matter-bounded and radiation-bounded compo- 
nents. 

The photoionization models do not provide a good fit to 
the observed optical spectrum of cloud 1. Given the prior 
evidence th at shocks contribute to the io nization structure of 
the XNC (Ford et al."1985'; 'Cecil"1988'), we compared the 
emission-line fluxes of cloud 1 with MAPPINGS II shock 
model grids dPopita & Sutherland 1995, 1996). The MAP- 
PINGS II shockH-precursor models with velocities in the range 
Vj = 400-450 kms"', magnetic parameters B/n^/^ = 2- 
4 /iG cvtr'l^, and scaled nitrogen abundances of ^^-^3.5 provide 
a reasonable fit to the observed emissi on hnes in this cloud. 
Later Chandra X-ray observations by iTerashirria & Wilso J 
d2001h further support the conclusion that shocks are the dom- 
inant source of ionization at the location of cloud 1 . 

8. CONCLUSIONS 

The optical and radio morphologies in the NLR of M5 1 ex- 
hibit a striking correspondence. The conical, limb-brightened 
optical morphology south of the nucleus is suggestive of evac- 
uation by the weak southern radio jet. The near coincidence 
of the bright optical emission clouds with the position angle of 
the radio jet and the proximity of the optical and radio knots 
also suggest that the radio jet may influence the kinematic and 
ionization structure of the optical emission. 

We find that, in general, photoionization is the dominant 
excitation mechanism for the brightest optical emission-line 
clouds. However, shock-nprecursor models provide the best 
fit to the optical emission-line cloud that lies near, but west 
of, the entry point of the jet into the diffuse radio lobe. Thus, 
shocks are an important excitation mechanism for at least 
some regions in the NLR of M51. Furthermore, we find 
some evidence for the presence of weak, higher velocity dis- 
persion clouds. Several [Olll] clouds were not well fit by 
a single Gaussian and required a second, broader Gaussian 
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emission component to obtain a good fit. The signal-to-noise 
and general spectral resolution of our data were insufficient 
to permit a more detailed analysis of this diffuse emission. 
Higher signal-to-noise, higher spectral resolution STIS data 
(Ferruit et al., in preparation) that map the optical emission 
with four slit positions will help determine if shocks play a 
more pervasive role in the NLR of this Seyfert galaxy. Ad- 
ditionally, higher signal-to-noise and higher spatial resolution 
radio maps would be beneficial in determining the importance 
of weak radio jet structures, observed in several Seyfert galax- 
ies, to the general AGN paradigm. 
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